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1 Introduction 

After some erratic appearances in the eighties, [1],[2] and [3], gravitational Chern-Simons 
(CS) action terms are met more and more frequently in the literature. This refers not 
only to a reappreciation of the earlier proposals (see, for instance, [4], [5]) but, also, to 
more and more frequent appearances of CS terms as part of the gravity action in various 
dimensions. Low-energy effective actions coming from superstring theories provide many 
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such examples. These terms are particularly interesting in relation to black holes. One of 
the problems they raise is how their addition modifies black hole solutions and associated 
charges. Another important problem is how to compute the black hole entropy in their 
presence. This is precisely what we wish to address in this paper. 

While the 2+1 dimensional case can count by now on a considerable number of analyses 
[6-10], little is known about higher dimensional gravitational Chern-Simons theories. In 
this paper we would like to start searching systematic answers to the questions raised by 
the presence of these terms in higher dimensions. We will extend our analysis also to 
reducible CS action terms containing non gravitational components. 

To be less generic, let as introduce a few definitions and basic properties. In this paper 
we wish to investigate some properties of gravitational actions extended with Chern-Simons 
terms, 

£ = £ cov + £cs (1.1) 

By C cov we denote some generic manifestly diffeomorphism covariant gravitational La- 
grangian in D dimensions, while £cs contain Chern-Simons terms and, therefore, is not 
manifestly covariant. 

Formally, a purely gravitational Chern-Simons term in D = 2n — 1 dimensions is 
obtained by going to (D + 1) dimensions, where 

P n (R,...,R) =dYg , (1.2) 

holds. Here P n denotes an invariant symmetric polynomial of a given Lie algebra (see 
Appendix A). We recall that for classical groups invariant polynomials can be defined via 
symmetric traces (str) of the Lie algebra elements. The P n 's we have just introduced for 
purely gravitational CS terms are the invariant polynomials with respect to the Lie algebra 
of the SO(2n — 1) group. Eq.(1.2) can be explicitly integrated (transgression) [11], 

Tg(r) = n/ dt P n (T,R t ,...,R t ) (1.3) 
J o 

The expression (1.3) will be taken to define an action term in D = 2n — 1 dimensions, i.e., 

Lpcs = Tg (1.4) 

where subscript pCS denotes contribution of the pure gravitational CS term. We will often 
drop the superscript W whenever no confusion is possible. 

It is important to recall that for n = 2k— 1, that is, in D = 4k— 3 dimensions, irreducible 
invariant symmetric polynomial vanish identically, Pik-\ = 0. So, purely gravitational 
irreducible CS terms can appear only in D = 4k — 1 dimensions. Generically P n will not 
be irreducible. 

We will also consider action terms where a gravitational CS term is multiplied by an 
invariant polynomial made of one or more gauge "curvatures", so as to fill up a D-form. 
Typically such reducible action terms will have the form 

Li, mix = Tg ) (r)P fc (F) or L 2jmix = P m (R)T[g(A) (1.5) 
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where F represents the generic curvature of the gauge connection A and is an invariant 
polynomial of order k, such that D = 2m + 2k — 1. A can be either a non-Abelian 
or an Abelian gauge connection, or even a RR field. We will refer to terms like (1.5) 
as mixed Lagrangians. Let us also notice that the two Lagrangian terms (1.5) descend 
via transgression from the same term P m (R)Pfc(F) in (D + 1) dimensions. We can have 
in general many of these fields simultaneously but, for simplicity, we will consider only 
one connection. The extension to several connections and curvatures is straightforward. 
Generally, one can have* 



In our subsequent calculations the various independent pieces of the Lagrangian yield 
independent contributions to the intermediate and final results. As a consequence they 
can be dealt with separately (this is generically not true for the equations of motion, 
but, as we will see, the latter need not be explicitly calculated). Therefore to make the 
presentation more clear, we will mostly use pure gravitational CS term L p cs to demonstrate 
the procedures. Calculations for the mixed CS terms usually proceed in a similar fashion. 
The total contribution to the results are obtained by summing the contributions of all the 
separate terms appearing in the total Lagrangian. 

A few comments on general CS terms are perhaps useful. They were introduced in 
the mathematical literature by S.S. Chern and J. Simons in 1974 (see [15]) to represent 
characteristic classes of the total space of a principal bundle. Gravitational CS terms have 
been introduced in theoretical physics mostly to describe topological theories. One should 
not forget that in Lcs(T), as well as in (1.5), there is no explicit dependence on the metric 
- the metric appears only through T. It is worth recalling that, while general gauge CS 
terms are not globally defined on the base space of a principal bundle (because a gauge 
connection A is generally not basic), the gravitational CS terms are globally defined in 
space-time due to the existence of a natural lift in the frame bundle, see for instance [16]. 

After these somewhat lengthy preliminaries let us come to the subject of the present 
paper, where we are interested in studying some consequences of introducing in a theory 
additional CS terms, in particular in the modifications induced in asymptotic or bound- 
ary charges, with emphasis on black hole entropy. For manifestly covariant theories with 
covariant Lagrangians, there is a powerful method based on the covariant phase space for- 
malism, which was adopted by Wald [17] and led to an elegant general formula for black 

*It can be shown that Lagrangians in which CS terms appear inside the definition of gauge field strength, 
like H = dB + Lcs in heterotic string theory, can be put in a classically equivalent form in which CS terms 
appear only through separate terms of the form (1.5), but not inside the definition of gauge field strengths. 
For derivation see, e.g., review [12], and for applications in calculations in all orders in a' see [13, 14]. 





where 
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hole entropy. In [18] the author outlined an extension of this elegant formalism to theo- 
ries which are not manifestly covariant (Lagrangians are not covariant) but have covariant 
equations of motion, such as (1.1), and applied it to the calculation of the black hole en- 
tropy in some instances when a CS term is present. In this paper we wish to carry this 
program to completion, so as to be able to apply it in general. More precisely we will give 
a more detailed and elaborate analysis of the general method outlined in [18] and apply it 
to general theories with CS terms in the Lagrangian (1.1). In the process we will introduce 
some necessary modifications with respect to [18] in order to guarantee consistency. We 
will compute in detail the CS induced modifications to Wald entropy formula, and find out 
that the final formula is very similar to Wald's original formula valid for covariant actions, 
with a significant modification. We will analyze the subtleties connected to covariance and 
to the use of different types of coordinate systems. Although the modified Wald formula 
for entropy looks non-covariant, we will show that it can be 'covariantized'. 

The paper is organized as follows. In the next section we introduce the main formulas 
and derive the generalized Cotton tensor and the symplectic potential. In section 3 we 
generalize the covariant phase space formalism to the case of actions containing CS terms. 
In section 4 we derive the entropy formula and discuss its covariance. Section 5 is for the 
conclusions. To render the paper more readable many of the calculations are postponed to 
dedicated appendices. 

2 Equations of motion 

Let us consider first the pure gravitational CS term in the Lagrangian (1.4)-(1.3). In order 
to compute contribution to the equations of motion we calculate the (generic) variation 
with respect to T. Using 



STt 



tST 



6K t = d5T t + [T t ,5T t ] = -D t 5T t 




we have 



<5L pCS (r) 



n [ dt(P n (8T,-R t ,...,Rt) + (n 





Jo 



Using (A. 8) this can be put in the following form, 



<5L pCS (r)=n / dt(p n (5T,R t ,...,H t )-(n-l)dP n (T,5T t ,K t ,...,K t ) 



■Jo 



+ 





It is easy to see that 



(n - l)P„(D t r, ST t , Rt, • • • , Rt) = (n - l)P n (d t Kt, 6T t ,K t , . . . , R t ) 



37-fn(£Tt> Rtj • • • j Rt) — -fn(5r, Rj, . . . , R() 
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This allows us to directly integrate one of the terms, obtaining finally 

5L pC s(r) = n (p n {ST, R, . . . , R) - (n - 1) d£ dt P n (T, 6T t ,H t , R*)) (2.4) 

By adopting a compact notation we can write this as 

<5L pCS (r) = nP^ffR"- 1 ) + de nc (T,ST) (2.5) 

where 

nc = -n(n-l) [ dt P n (r,5r t ,R™- 2 ) (2.6) 
J o 

It is immediate to generalize (2.4, 2.5) to mixed Lagrangians, 

5Li, mfa = mP m (5T, R^P^F) + *P m (R)P fe (<5A, F^" 1 ) - (2.7) 
-d(m{m-l) J dtP m (T,5T t ,K?~ 2 )P k {F) + kr { ^{T)P k (5A,¥ k 



5U, mix = mP m (,5r,R m - 1 )P fc (F) + feP m (R)P fe (5A,F fe - 1 ) + (2 
+d (mP m (5r,R m - 1 )Tg(A) - k(k - l)P m (R) J dtP k (A, 5A t , F f fe ~ 2 



Let us concentrate first on (2.4). We introduce the 2(n — l)-form K = R n 1 ; its 
components are 

jy-afS _ (2ra — 2)! j o CT n-2^ (9 n\ 

This means that we have 

5L pCS = n(2n - 1)K V, 2 ... M2 „_ 2 ^ • • • dx^ + dO nc (2.10) 

Using the Levi-Civita tensor the previous formula can be rewritten as 
n( — 1 ) s 

OJjpCS - 7^73^7 6 A/i. lA12 -M2n-2 ( ' i QM2n-l V 5 « » + « U l^ 11 ] 



(2n - 2) 

where s denotes the metric signature. Now, using 

^ = \ 9 Xf) (V a <^ + V M % Q - 5g a ^) (2.12) 
and Bianchi identity we obtain 

«5L pCS = C^Sgfr e + dO cov + d@ nc (2.13) 
where e = y/—g d 2n ~ 1 x is the volume D-form, the components of COV are 



C °Vm2....2„-2 = (-l) s n(*KY^5g^e c 



/JlA t 2-"A t 2n-2 V > \ I VPt 1 &/J.lH2---fl-2n-2 

i K al3 

\a\ll2—H2n-2] 



to 2 6g a . a ,K*.., (evenn) (2.14) 
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and 



= -(-l) s nV a (*]f) a(W (2.15) 

where *K is the Hodge dual of K: 

y*^> ~(2n-2)! mv2-H2n-2 € t z - w ) 

Obviously is the contribution to the equations of motion of the L p cs term. It is 
a symmetric, traceless and divergence free (pseudo)tensor, investigated in some detail in 
[19] where it was shown that it can be viewed as a generalization of the Cotton tensor to 
D = 2k + 1, k > 1 dimensions [19, 20]. 

As can be read off from (2.13), the total symplectic potential & is the sum of the two 
terms COV (given by (2.14)) and nc (given by (2.6)). We refer to G cov and nc as the 
covariant and the noncovariant contribution to the symplectic potential, respectively. 

3 Covariant phase space formalism for noncovariant Lagrangians 

The covariant phase space formalism is a powerful tool to calculate various charges con- 
nected with asymptotic symmetries. For example, it was used in [17, 21, 22] for an elegant 
derivation of the first law of black hole thermodynamics in theories with manifestly diffeo- 
morphism covariant Lagrangian descriptions. In [18] the formal extension to theories with 
Lagrangian description which is not manifestly diff-covariant was outlined. In this section 
we wish to apply this construction to theories with general Chern-Simons terms. In the 
process we calculate all the necessary ingredients that will be used in the next section to 
obtain a general formula for black hole entropy in such theories, extending in this way the 
results of [17, 18]. 

In what follows, as usual, we understand that the Lagrangian is a D-form. Also we 
denote by the symbol a (j>" all the dynamical fields, that is the spacetime metric g a b as 
well as any matter field. The derivation of the first law in [21] assumes a diff-covariant 
Lagrangian, for which the following holds, 

<%L cov ((/)) = £^L COV (0) 

where 6e is the variation induced by the diffeomorphism generated by the vector field £ 
with components £ a , while £^ denotes the Lie derivative with respect to £. This condition 
is not satisfied if there are Chern-Simons terms in the Lagrangian; instead, we have 

<%L(</>) = £{L(<f>) + dSz (3.1) 

for some (D — l)-form 3^. The form of (3.1) expresses the fact that the action is still 
diff- invariant (if the space-time is closed), though the Lagrangian is not diff-covariant. 
It is understood that £^ acts on non-tensor quantities (such as T) as if their indices 
were tensorial ones. As already pointed out above, in [18] the problem of extending the 
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covariant phase space procedure to Chern-Simons type Lagrangians was first dealt with. 
We will follow in part this recipe and introduce the necessary modifications. 
The first-order variation of the Lagrangian L takes the general form 

SL(4>) = E<J0 + d0(^,(50) (3.2) 

where the D-form E represents the equation of motion ^ (which is manifestly diff-covariant) 
and the 'surface term', (D — l)-form 0, is known as symplectic potential, as we have already 
pointed out. In Sec. 2 we have already found that CS terms in the Lagrangian introduce 
additional contribution to the symplectic potential which can be separated in covariant 
and noncovariant part 

= @nc + qcov (33) 

where in the case of pure gravitational CS Lagrangian term (1.4), nc and COV are given 
by (2.6) and (2.14), respectively. 

For the two mixed CS terms (1.5) we obtain 

&l c (T,5T,A,5A) = -m(m-l) f dt P n (T, 5T, R™~ 2 )P fc (F) - 

Jo 

- fcT< ro >(r)iV$A,F*- 1 ) (3.4) 

and 

&% C {T,5T, A,5A) = mP m (5r,R m - 1 )T( fc) (A) - 

- k{k - l)P m (R) / dtP k (A,5A t ,Ft 2 ): (3.5) 
Jo 

respectively. The covariant parts can be easily obtained, but, as we shall see, they are not 
essential for deriving the black hole entropy and so we shall not bother to write them down. 
Next, we define the non-covariant part of the diff transformation by 

<% = - £ i (3.6) 

Let us recall that for tensor- valued p-forms (like 1-form T and 2-form R) the Lie-derivative 
satisfies 

i? s = di € + i s d+[A, ] , A b a = d a £ b . (3.7) 

In particular, 

= £^T + 5^T , <5 c r = dk (3.8) 
This allows us to define the Eg term via the relation 

<5 € L(r) = d%(r) (3.9) 

' Summation over the dynamical fields and contraction of their tensor indices with corresponding dual 
tensor indices of E are understood. 
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For the pure gravitational Chern-Simons Lagrangian (1.4) we get 

H^ n) (r) = n(n- 1) / dt(t-l)P n (dA,r,R t "- 2 ) (3.10) 
Jo 

For the mixed Lagrangian case we have simply 

S u (r,A) = 3 { f\T)P k (F) and E 2 , e (I\A) = (3.11) 

One basic ingredient we need is the Noether charge Q^. This is defined via the Noether 
current (D — l)-form J^, 

J 5 = 5 - z f L - Eg (3.12) 

which is conserved on-shell (that is, for the field configurations that satisfy E = 0), 

dJ^ = dO^ — d{i£~L) — dS^ « <5gL — d(igL) — dS^ = £^L — d{i^L) = 0. 

where ~ denotes equations that are valid on-shell. Here we have also used (3.7), which for 
a Lagrangian L reduces to 

£^ = di£ + i£d (3.13) 

and the fact that dh = since it is D-form. This implies, according to [21, 23], that is 
exact on-shell, i.e. there is {D — 2)-form Qg, such that 

J 5 « dQ ( (3.14) 

We refer to this property as Wald lemma. 

To find Q^, as a local expression of the dynamical fields linear in £, one can use either 
the constructive method from [23], or more formal methods based on cohomology and the 
so-called variational complex. In all our cases Qg takes the form 

Qe = Qf } + Qf (3.15) 

where the superscript denotes the number of derivatives acting on £. The details of calcu- 
lations are presented in Appendix C. 

The results are as follows. For the purely gravitational case we obtain for 



/' 

JO 



QY> (V) = n{n - 1) / dtP n (A, T, R" ) (3.16) 



and for (in components) 



(Q 



fW"A*»-a =<n-\) / dtP n (^r,r t ,Rr 2 ) M1 ... M2n _ 3 

J o 

-nH b K ab a ^ 2n _ 3 + ^(2n - l)^[a^£,.. /i2Tl _ 3 ] (3-17) 
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For the first mixed CS term we have 



Q^(r,A) = m(m-l) I' dtP m (A,T,R?- 2 )P k (F) (3.18) 

J o 



and for the second: 



Q«(I\ A) = mP m (A,R m " 1 )Tg(A) (3.19) 



As we will explain later, only the part is relevant to the construction of the black hole 
entropy. 

Proceeding further we can calculate 55^ L in two ways, 

<%5L » 5 c d0 = d<%0 = d(^0) + d<%0 = di^d® + d<%0 
By comparison we see that 

d<%0 » 5dS { = d5E { (3.20) 
having used [5, d] = in the last equality. This allows us to introduce a (D — 2)-form 



via 



<5 C - 63$ « cffi 5 (3.21) 

With the usual methods, for the purely gravitational CS Lagrangians, one can compute 
(see Appendix B.2) 

S f (A,r) = -n(n-l)(n-2) / (it t(t - l)P„(dA, T, (5T, R^ 3 ) (3.22) 

We remark that, actually, in this case we have an exact relation: 

<%0 - <53 f = effif (3.23) 
In a similar way, for the first mixed Lagrangian we have the exact relation 

<%©r (r, A) - 5H 1>e (r, A) = ds 1>e (r, A) 

where 

S u (r,A) = -m(m-l)(m-2) / dt(i 2 — t)P m (dA, T, 5T, R™~ 3 )Pfc(F) 

JO 

-fcm(m-l) / dt(t-l)P m (dA,r,R t m - 2 )P fc (5A,F fc ~ 1 ) (3.24) 
■/ o 

Finally, it is easy to prove that S 2 ,g = 0. 

A central element in covariant phase space approach is the symplectic current (D — 1)- 
form u defined by 

u>(4>, 6 1( />, 5 2 0) = 5i0(</>, <5 2 0) - <5 2 0(0, Si<f>) (3.25) 
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Let C be a Cauchy surface with orientation given by 

where n a is the future pointing normal to C and e is the positively oriented spacetime 
volume form. We define the presymplectic form $7, a 2-form in the space of fields and 
0-form in spacetime, by integrating the (D — l)-form u over a Cauchy surface C, 

Q(<P, 5 2 (f>) = J u>(4>, S 1( j), 8 2 4>) (3.26) 

The Hamilton equations of motion generated by vector field £ are then 

SH[^]=n((j),84>,6 $ (j)) (3.27) 
where H[£] is Hamiltonian for £. From (3.25) and (3.12) one obtains (see [18]) 

<5J£ pa 6<j>, Stf) + c%0 + S e ) (3.28) 
which, by using (3.26-3.28), leads to 

« / (<5Q S - z 5 e - £ e ) (3.29) 

We see that the existence of H[£] requires either that the last two terms on the right hand 
side of Eq. (3.29) give vanishing contribution, or that there is a (D — 2)-form satisfying 

SCz = ^0 + E ? (3.30) 

If the latter were the case one could introduce = — and integrate Eq. (3.29) to 
obtain 

H[£] ~ / (3.31) 
JdC 

However, even for the pure gravitational Chern-Simons Lagrangian of the form (1.4) it 
seems to be generically very difficult to construct a Cg satisfying (3.30) off-shell. We 
therefore turn to the first alternative, which requires studying the problem on a case by 
case basis [24]. 

In our application to black hole entropy our strategy will be the following one: first 
perform the calculations in Kruskal-type coordinates, then show that i^Q + vanishes, 
and finally "covariantize" the results. 

4 Black hole entropy in theories with Chern-Simons terms 
4.1 BH Thermodynamics in Kruskal coordinates 

Wald and others ([21], [22], [17]) have shown how to compute the black hole entropy using a 
covariant phase space formalism. They have considered a general classical theory of gravity 
in D dimensions, arising from a diffeomorphism-invariant Lagrangian. Assuming that the 
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2tt I d D - 2 x^-g— e ab e cd (4.1) 

E 



theory admits a stationary black hole solution with a bifurcate Killing horizon and that 
the canonical mass and angular momentum of the solution are well defined at infinity, the 
first law of black hole mechanics always holds for perturbations of nearly stationary black 
hole solutions. Using this, it was shown ([17]) that for a given Lagrangian C, the entropy 
of a black hole can be computed by an elegant formula 

5C_ 

SRabcd 

where the integration is over the horizon cross section £ with binormal e ab normalized as 
£a&e ab = — 2. We now want to find a generalization of Wald entropy formula (4.1) to the 
theories with Chern-Simons terms of the most general form. 

Following [17, 18] we assume £ to be a Killing vector field, and that dynamical fields 
are symmetric, i.e., 5^<p = 0. Then (3.29) gives 

(5Q 5 - z ? - £ e ) m 6H[£] » (4.2) 

dC 

For an isolated black hole the boundary of the Cauchy surface dC consists of two discon- 
nected parts, one is an intersection of C with the horizon, which we denote by S, and the 
other is the asymptotic infinity (denoted by oo). So (4.2) becomes 

/ (5Q 5 - ^0 - S f ) « / (6Qz - i£® - S 5 ) (4.3) 

JE Joo 

As in [17] we assume the black hole has Killing horizon, generated by £, and that the surface 
gravity k is constant on the horizon. * One consequence is that such spacetime admits an 
extension in which the horizon has a bifurcation surface B [26], on which 

e| B = (4.4) 

This property makes it extremely useful to choose the Cauchy surface such that £ = B. 
From now on we are assuming this. 

Now we want to show that (4.3) yields the first law of thermodynamics, namely 

T5S = 5U - ... . (4.5) 

The left hand side of (4.3) needs a more careful analysis, especially as the individual 
terms are (at least apparently) not covariant. In "regular" coordinate systems (which are 
well-defined at the bifurcation surface B), property (4.4) implies that the second term in 
(4.3) vanishes 

it@\ B = 0, (4.6) 

From the same argument it follows that we can drop the part of Qg. Also, in coordinate 
systems regular at B, we can write A in the following covariant way 

^\\B = ^bt a \ B = Ke\ (4.7) 



* It was shown in [25] that k is constant on Killing horizon for all stationary-axisymmetric black holes 
with the "t — <fi" reflection isometry. 
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s 



where e is the binormal of B. 

We still have to deal with the second term in (4.3). As we do not know how to 
generically integrate Xf , we turn to a different strategy. From the results of the previous 
section we saw that Xg contains a dA factor. Consequently it appears beneficial to specialize 
for the moment to a particular "Kruskal-type" coordinate system (details are in Appendix 
F) in which 

dA = . (4.8) 

This makes the Xg term trivially vanish^. 

From the form of , obtained in the previous section, and (4.7), it follows that the 
left hand side of (4.3) can be written as 

/ Qt = ^5S hh = T6S hh (4.9) 

As in manifestly covariant theories, the black hole temperature T can be identified with 
k/2it. The black hole entropy S"bh has two contributions, 

S hh = S cov + S cs (4.10) 

where S cov comes from the manifestly covariant part of the Lagrangian L cov , and Scs from 
the noncovariant part of the Lagrangian Lcs- S cov is obtained by putting C cov in Wald's 
formula (4.1), while Scs can be obtained by using calculated in Section 3. For example, 
from (3.16) follows that the contribution from the pure gravitational Chern-Simons term 
(1.4) to the black hole entropy formula is 

S CS = 2vrn(n - 1) f dt [ P n (e, T, R™~ 2 ) (4.11) 
Jo Jb 

It is interesting that C-S contribution to the entropy formula (4.11), and analogous 
relations for mixed cases, can also be put in the form analogous to the Wald formula (4.1). 
To see this let us write Wald entropy formula (4.1) in the equivalent way^ 

Sw = 2. [ e\ ( (4.12) 



B \ SRa bvv 



VP1---PD-2 



§ What is at stake in relation to the term (which is anyhow absent in 3D) is the following. If is 
not integrable, i.e. there does not exists any local Cj such that = SC^, then a charge can still be 
formally obtained by integrating along a path in the space of solutions, starting from some conventional 
one. However in this way the result turns out to be path dependent. Path dependence of this result and 

non-integrability are two equivalent statements. In this situation the entropy would not be uniquely 
defined. It follows that we are obliged to assume Sj-integrability. However, after assuming it, we still do 
not know the explicit form of Cf, which is very hard to construct. The procedure described in this section 
teaches us how to avoid such explicit construction. 

^This follows from 

£d L = e A Zd-2 

where eo is the volume D-form, <Ln-i is the induced volume (D — 2)-form on B, and e is 2- form binormal 
to B. 
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For C-S contribution to the entropy we obtain 



fiVpl—pD-2 



'dn t % J B "\8Rt- blw 



and we see that the only difference with (4.12) is that one has to consider the variation 
with respect to Rt instead of R. Here it is understood that variation with respect to R< 
acts inside the t integral present in Lagrangian. The Lagrangians used here are of the 
form LgQg from (1.6) which means that there will always be exactly one T factor, and 
consequently exactly one t integral. 

As we show in Appendix F (see e.g. relation (F.14)) using Kruskal-type coordinates 
enables us to substitute Rf with iR in (4.11). We then obtain 

5 CS = 2vrn f P n (e,T,R n ~ 2 ) (4.14) 
JB 

In the mixed case 1 (1.5) one gets 



Si.mixcs = 2vrm I P m (e, T, K m ~ 2 )P k (A) (4.15) 
JB 

2vrm / P m (e,R m - 1 )Tg(A) (4.16) 

JB 



IB 

<S2,mixCS 



IB 

In the above integrals it is of course understood that T, R and A are replaced by their 
pull-backs by the immersion of B in space-time. The pull-backs of Y and R have in general 
components in all of the Lie algebra of SO(l, D — 1). However the use of Kruskal-type 
coordinates, adapted to the black hole geometry, entails that the only contribution to 
such terms comes when the tensor indices are in the 2-dimensional normal bundle of the 
[D — 2)-dimensional bifurcation surface B. In other words 

str(erR 2fc ) = 2T N ~R$ , (4.17) 

where 

r w = ^tr(er) , Rjv = ^tr(eR) (4.18) 

Later on we shall use (4.17) to analyse the issue of coordinate independence of black hole 
entropy in theories with CS terms in the action". 

To derive (4.17) we note from (F.4) and (4.7) that the components of the binormal in 
Kruskal coordinates are 

e u u = l, e V v = -1 (4.19) 
then by using (F.21)-(F.22), it follows 

str(erR n " 2 ) = Y u u{Yi u u) n - 2 - T v \ (R v ' v ) n ~ 2 



"Formula (4.16) is of the covariant type entropy formulas, since A is simply a one- form from the diff 
point of view. Therefore from now on, when talking about covariantization, we will refer only to formulas 
(4.14,4.15). 
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On the bifurcation surface we have R u u = -B/ v = H N . Since n is even it follows 

str( e rR"- 2 ) = (r u u-T v v ){iL N ) n ~ 2 

= 2T N R n N 2 (4.20) 

where T and R denote the pull-backs on B. The 'miracle' of Kruskal-type coordinates 
is that (locally) they kill all the components of T, except those in the normal bundle 
directions. 

Using (4.20) the formula for the entropy becomes 

5 CS = 47rn I T N R™" 2 (4.21) 
Jb 

4.2 Covariance of black hole entropy 

In Section 4.1 we have obtained a formula for the black hole entropy by using a special type 
of coordinate system - Kruskal-type coordinates. This raises immediately the question of 
the validity of the same formula in a generic coordinate system. Beside, in practical cal- 
culations one rarely uses Kruskal-type coordinates, and "Schwarzschild-type" coordinates 
(which, on the other hand, are not even regular at B) are typically the most efficient. In 
reality it is possible to derive the same formula (4.21) in a purely geometric and coordi- 
nate independent way: see Appendix E. The geometry we have to do with in this section 
is determined by the presence of a submanifold B, which naturally breaks the group of 
diffeomorphisms into the direct product of Diffjg and the diffeomorphisms that leave B 
pointwise unchanged. The latter can be interpreted as gauge transformations of the B 
normal bundle, NB. For our problem the relevant invariance to consider is the invariance 
of the entropy formula under this product group**. 

Although this problem can be dealt with in more generality (see Appendix E), it is 
interesting to examine the question of the black hole entropy formula(s)' invariance from 
the point we left it in the previous subsection, i.e. the more popular point of view of 
Kruskal coordinates. Our purpose in this subsection is to 'covariantize' (4.21), that is to 
find a covariant formula for entropy that coincides with (4.21) in Kruskal coordinates. This 
covariant formula will turn out to coincide with that found in Appendix E. 

In the previous subsection we have obtained that the contribution to the black hole 
entropy due to the presence of CS terms in the action, calculated in Kruskal coordinates, 
will be in the most general case a sum of terms of the form 

2vr(2A:i+2) / str( e rR 2fcl )str(R 2fc2 ) • • • str(F mi ) • • • , 
Jb 

A ^2 k i + 2 ^2 m j = D ~ 3 ( 4 - 22 ) 

* 3 

We have seen that this can be reduced to^ 

2vr(2A;i + 2) [ T N R 2 v fcl str(R 2fc2 ) • • • str(F mi ) • • • (4.23) 
Jb 

"One can consider also more general diffeomorphisms, see for instance [27], but we will ignore such 
kinds of subtleties here. 

^Once again we recall that by T and R in (4.22) we mean the pull-back to B of the corresponding forms 
in space time. 
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Now, in Kruskal coordinates we have 

qlY Nv = q v tlU PV v v p = ^ (4.24) 

where and are null vector fields (normalized to u^v^ = — 1) normal to B defined in 
(F.7) in Appendix F, and q„ is a tangential projector on B (first fundamental tensor of B) 
defined in Appendix D. The second equality is simply the definition of u, which is nothing 
but the connection of the normal bundle, see (D.10). 

Next we express also Rat in terms of u. Using (D.15) and (D.12) we get 

qP qZ (Rat) pct 1 fi M „ = qP q°(du) pa (4.25) 

where fi is the normal curvature 2-form. 

What is important here is that: (a) by definition, u is a 1-form, fl is a 2-form on B 
(b) inside the integral (4.23) one can substitute T^p with q^Y^u-, and also (Rat) mi/ with 
q P l q"{R.N)pcT-, because integration over B implies tangential projection by q^. In particular 
we remark that u behaves as a 1-form under the diffeomorphisms of B and d is the exterior 
derivative on B. So the integral (4.23), with V replaced by u, is well defined and covariant 
with respect to the diffeomorphisms of B. In summary, we should replace (4.23) with 

27r(2fci + 2) / ^fl 2fcl str(R 2fc2 )---str(F mi )--- (4.26) 
Jb 

and use this form of the integral to prove its invariance under normal bundle gauge trans- 
formations. This is simple. The structure group of the normal bundle is SO(l,l). The 
corresponding gauge transformations are given by rescalings —> fv^, — > u^/f (see 
Appendix D), where f{x) is an arbitrary positive function. From (4.24) it follows that u 
transforms as 

u-^u + dlnf (4.27) 

From this the transformation rule 

u{du) 2k u(du) 2k + d(. . .) (4.28) 

follows. Finally, we see that, as a consequence of (1.2), the formula for entropy (4.26) is 
invariant under the gauge transformation (4.27). 

In conclusion, we have shown that the expressions for the CS black hole entropy terms 
(4.26) are invariant. 

4.3 BH entropy in Schwarzschild-type coordinates 

In the previous subsection it was shown how the contribution to the black hole entropy 
formula from the CS Lagrangian terms can be written in the manifestly diff-invariant 
form. Now we turn to the question of practical calculations. In practice, calculations are 
typically performed in coordinate systems in which symmetries are maximally exploited, 
which typically means that a maximal set of linearly independent Killing vector fields are 
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used as generators for coordinates. We call such coordinate systems "Schwarzschild-type" 
(because the coordinates used in the Schwarzschild solution are one such example). The 
bonus of these coordinate systems is that if there are iV Killing vectors, then the components 
of the metric depend only on (D — N) coordinates, which substantially simplifies the 
calculations (in particular, to solve the equations of motion). In our calculations, the 
obvious simplification comes from the property that in Schwarzschild-type coordinates, by 
construction, the components of all Killing vectors are constant everywhere, and this is true 
in particular for the horizon generating Killing vector £ a . This means that A a {, = db£, a = 0, 
and so also dA = 0, everywhere. This implies T,^ = 0, like in Kruskal-type coordinates. 

However, there is a price to pay. Schwarzschild-type coordinates are by construction 
singular on the bifurcation surface B (in fact, on the whole horizon %). The constancy of 
the components £ a , due to (4.4), implies that the metric tensor components, calculated in 
Schwarzschild-type coordinates, must be singular on B. This means, for example, that (4.7) 
and (4.6) are not valid in these coordinates. Therefore one has to be careful when using 
Schwarzschild-type coordinates to compute quantities at the horizon, and the entropy in 
Wald formalism is one such quantity. In manifestly covariant theories this singular behavior 
does not pose a problem - physical quantities typically can be put in diff-invariant form 
which allows one to perform calculations slightly outside the horizon and then take the 
horizon limit which is well-defined for diff-scalars. When CS terms are present, manifest 
diff-covariance is broken and, as a consequence, one has to make sure how and when one can 
calculate with Schwarzschild-type coordinates. One of the benefits of the covariantization 
from the previous subsection is that we can safely use Schwarzschild-type coordinates in 
calculating the entropy terms in the form (4.26). 

What we want to show here is that in Schwarzschild-type coordinates the CS contri- 
butions (4.26) to the entropy can be calculated directly from 



We shall now prove this explicitly for the case of stationary black holes. Schwarzschild-type 
coordinates, and their connection to Kruskal-type coordinates, are presented in Appendix 
G. It was shown in [26] that "natural" (fixed-time) Cauchy surfaces Ct (defined by t = 
const) are ending on the bifurcation surface B, which is convenient for analysis of the first 
law of thermodynamics, which follows from (4.2)-(4.3). As we show in detail in Appendix 



where in subscript we denote the coordinate system in which the calculation is performed. 
As all other factors in (4.29) are manifestly covariant, from (4.30) it follows that it does 
not matter whether we calculate (4.29) in Schwarzschild-type or Kruskal-type coordinates. 
This completes the proof. 




(4.29) 



G 



(FaO Schwarzschild _ (rV)Kruskal 



(4.30) 
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5 Conclusion 



In this paper we have shown that the addition of CS terms to a gravitational action 
can be dealt in a clear, if not straightforward, way. We have extended the covariant 
phase space formalism to such noncovariant situations and found a general formula for 
the black hole entropy. This formula looks very similar to the covariant Wald formula, 
but with a significant difference that seems to render it non-covariant. By exploiting 
the geometry of the bifurcation surface we have however been able, by means of Kruskal 
coordinates, to reduce the entropy formula to a form that can be 'covariantized'. In other 
words we have found a covariant formula that reduces, in Kruskal coordinates, to the 
latter. The final 'covariantized' formula is covariant with respect to the diffeomorphisms 
of the bifurcation surface and to the gauge transformations of the normal bundle. We have 
outlined a derivation of this formula from the initial seemingly non-covariant form in a 
coordinate independent way, via the geometry of the relevant orthogonal bundles. 
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Appendix 

A Mathematical preliminaries 

We denote q- valued (q being a Lie algebra) p- forms as follows: 

A = A a T a = L A a ^,„^T a dx 1 " 1 A--- Adx>* E q®K p {^) (A.l) 

Commutators between such objects are defined in the following manner 

[A, B] = A a B b [T a , T b ] = AB - (-1)™ BA (A.2) 
where A a is p-form and B a is g-form. In the case of the Christoffel 1-form this means that 

[T, T]\ = 2(T 2 ) a b = 2T% Tl dx» A dx v (A.3) 
Throughout the paper we are using the standard notation, 

r = dv + r 2 , r t = tr , Rt = dv t + r 2 = m + (t 2 - t)v 2 (A.4) 
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Covariant exterior derivatives are 

B = d+[T, .] , D t = d+[T t , .} (A.5) 

and Bianchi identity reads DR = and D^R^ = 0. 
We often use the identity 

d t TLt = R + (2t - l)r 2 = dT + 2tT 2 = dT + [T t , T] = D t r 

and the variation of the Riemann 2-form, 

6Kt = 5dT t + 5T t T t + T t ST t = d5T t + [T t , 5T t ] = B t 5T t 

The noncovariant part of the diff-variation 5^ is defined by 

<5 5 r = dA 

and implies 

<%r t = tdA, <%R = , 5 e Rt = (t - l)(r t , dA] = (t - l)D t dA 



Let Mfc(C) be a set of complex A: x k matrices, and P n : <S> n Mjt(C) — >■ C n-linear functions 
(linear in each of its arguments), which are 

a) symmetric, 

-fn(-^-i) • • • i ^4ij . . . , Aj, . . . , j4 n ) = P n (Ai, . . . , j4j, . . . , Aj, . . . , vl n ) (A. 6) 

for all A r G M fc (C), 1 < r < n, and 

b) adjoint invariant with respect to the gauge group G, 

P n {g~ 1 A l g,...,g- 1 A n g) = P n (A l ,...,A n ) (A.7) 
for all so g £ G and A r E g = Lie G. 

A map P n that satisfies these properties is called symmetric invariant polynomial. An 
example of such polynomial is symmetrized trace, 

Str (A X ,..., A n ) = — ^ tr (Ar(l) • • • Ar(n)) 

7T 

where 7r denotes the permutations of the indices {1, . . . , n}. 

When some of the arguments of the invariant polynomial happen to be equal, we often 
use the following abbreviations, 

P ri (A,B"" 1 ) = P n (A,B,...,B) 
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P n (A, B, C n ~ 2 ) = P n (A, B, C, . . . , C) 
etc. 

It is easy to extend the domain of invariant polynomials from q to g-valued p-forms on 
manifold by 

P n {A 1 rn,...,A n T) n ) = P n {A 1 ,...,A n )r] 1 A--- Ar) n , rji e A Pi (M) 

We often use the following identity (proof can be found e.g. in [28]), 

n 

dP n (A im , A nVn ) = £(_i)w+~+«-i P n (A m , D(A iVi ), ...,A n ) (A.8) 

i=l 

where covariant derivative D is defined as 

D(A lVi ) = Aid m + [0, AiVi] = Aidw + 9*A\ [T a , T b ] dx» A m 
for some 3-valued 1-form 9. 

B Some exact formulae 

In this Appendix we derive some exact relations, which are used in the text or will be 
useful later on. Here exact means that they are true also off shell. The purpose of this 
appendix is also to show some of the techniques we use in deriving our formulas. 

B.l Noether current 

Let us start from 

J ? = 0^ + 0- - i e Lpcs - S 5 (B.l) 

where 

Q nc _ & nc( §T = ^ d + + DA j 

Using the methods from section 2 we have 

& nc (5T = DA) - S € = dq 14 - nP n (A, K^ 1 ) (B.2) 

where 

q u (r) = n(n - 1) f l dtP n (A, T, Rr 2 ) (B.3) 
J 

Now let us use 



i € R 4 = L € r t - £> t (**r t ) 



(B.4) 



where = i^d + dig. We have 

i ? L pCS = n f dtPni^T,^ 1 ) - n(n - 1) f dtP n (T, LgT t - D t (i 6 T t ),R« 
Jo Jo 

= nf <ftP n (i € r,R£ -1 )-n(n-l) [ dtP n (T, LgT t , R?" 2 ) 
Jo Jo 

f 1 f 1 d 

-n(n-l)d dtP n (i^T,T t ,B^~ 2 ) + n / dtt—P n (igT,K 
Jo Jo dt 

= nP n (^r,R n - 1 )+n(n-l) / dtP n (LgT, T t ,R 

Jo 

-n{n - l)d [ dtP n (igT, T t , R"'~ 2 ) 

JO 



r 1 ) 



On the other hand 



@ nc (ST = (igd + di € )r) = n(n - 1) / dtP n (L(:T, T t , R 

■/ o 



r 2 ) 



It follows that 



e nc (jr = (* e d + d»f)r) - ^l pCS = -nP n (^r, r"- 1 ) + d q2 , 5 



where 



q 2iC = n(n - 1) f dtP n {i^T, T t , R™~ 2 ) 

JO 



So, defining, 

q? (r) = q u (r) + q 2)4 (r) 



n(n — 1) 



/ dtP n (A, T, R™" 2 ) + / dtP n (igT, 
Jo Jo 



r*,Rr 2 ; 



we have 



J ? = 0«w _ nP n (i 5 r + A, R ) + dq c 



(B.5) 



(B.6) 



(B.7) 



(B.8) 



(B.9) 



The obstruction to the off-shell exactness of is Q^ ov — nP n (i^T + A,R n ). This 
obstruction is removed by putting the system on-shell (see below). This will allow us to 
define Qg. 

One can derive similar exact equations for the mixed cases. 
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B.2 Derivation of 



We give here the derivation of in the purely gravitational case. 

<%0 - <5H C (B.10) 
= (n(n - 1) J dtP n (T, ST t , R™~ 2 )^ - 5 (n(n - 1) J dtP n (dA, T, R™~ 2 )^ 

= - n (n-l)(jf diP n (dA,5r t ,R"- 2 ) + (n-2)(t 2 -t)P n (r,5r i ,[dA,r],R™- 3 ) 

-n(n - 1) ^ di (t - 1) (> n (dA, <5T, R™~ 2 ) + (n - 2)P n (dA, T, D t <5r t , R™" 3 )) 
= -n(n- l)(n- 2)d / dt (i - l)P„(dA, T, <$r t , R™~ 3 ) 

_ n(n _i)( n _2) dt(t-i)p n (dA,^,5r 4 ,Rr 3 ) 

-n(n-l) / dt(2t- l)P n (dA,<5r,R™~ 2 ) - n(n - 1) / di P„(dA, ST t , Bl 
Jo Jo 

= -n(n- l)(n- 2)d / dt (i - l)P„(dA, T, <$r t , R™~ 3 ) 

JO 

_ n(n _ i)( n _ 2 ) ^ 1 ^ _ i)p n (dA, «yr t , Rr 3 ) 

= -n(n- l)(n- 2)d / dt (i - l)P„(dA, T, <JT t , R: 
7o 



r 2 ) 



That is 

E f = -n(n-l)(n-2) / dt (t - l)P„(dA, T, <$r t , R™~ 3 ) (B.ll) 

JO 

The generalization to the mixed cases is straightforward. 

C Relations between No-ether current and charge 

In this appendix we present various ways of deriving the relation 

J 5 « dQ ? (C.l) 

from the on-shell closedness of J^, dJ^ ~ 0. In (C.l) both and are local expressions 
of the fields and their derivatives, and what is remarkable of the relation is of course the 
specification that it holds on shell. Without this specification the relation would simply be 
the statement of the local Poincare lemma. Modding out the equations of motion requires a 
more sophisticated version of the local Poincare lemma. This problem has been tackled and 
solved by Wald, [21], in a constructive way. But the same problem had been dealt with and 
solved in a more mathematical language, using the cohomology of the so-called variational 
bi-complex. In fact solving this problem is equivalent to finding the homotopy operator(s) 
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showing the exactness of such bi-complex. On the other hand Wald's constructive approach 
is tantamount to constructing the relevant homotopy operator (s). 

In this Appendix we follow a line in between. We use a constructive approach, by 
introducing a simple practical method to invert eq.(C.l). This method relies on a new bi- 
complex that can be embedded in the standard variational bi-complex, thus transforming 
it into an variational tri-complex. Its exactness completely justifies our heuristic approach. 

C.l Bicomplex formulation of Wald's Lemma 

In the cases of our interest, Jg is always of the form 

J f = dd a C b X\ + d a e + £ 6 A(°) 6 (C.2) 

where X, and A<°) are tensor- valued forms (which are local functionals of dynamical 
fields 4>). We know that on-shell (i.e., when <j) satisfy EOM) Jg is closed, so 

dJ c = =► J € = dQ(= (C.3) 

Our goal is to calculate Qg. 

We now show that Wald's constructive result [21] can be cast into an elegant coho- 
mological formulation. In this subsection we will simply ignore the specification that the 
relations we are dealing with are true on shell. This will be justified in the next subsection, 
by inserting the following construction into an exact variational bicomplex. 

Following Wald, we notice that in all formulae indices which are contracted with partial 
derivatives of £ appear only in symmetrized combinations, so we work only with tensor- 
valued p-forms of the form 

TT Ol-ttj , — TT. ,(oi-Oj) (0 41 

11 Hi— Hp b — n [ni-fip] b {^-^J 

which we denote in short Hp 1 , and in general call p-tensors. By definition, when j < or 
p < then Hp = 0. The exterior differential d changes p-tensor to p+1 -tensors. We define 
two additional operations on p-tensors. One is the "g- wedge" operation, g, defined by 

(gHt) n+ r~ ai+1 b = (p + l)g^H n+l] «*~^\ (C.5) 

(g® = <5° is the metric tensor) which transform p-tensors into P ^}-tensors. We also define a 
"^-contraction" operation, i g , by 

(i 9 ^) M ,.. Mp _ 1 ^'-^- 6 = i^(^)^ 1 ...^_ 1 --^-i fe (C.6) 

which transform ^,-tensors into P Zi-tensors. The anticommutator of g and i g satisfies 

{i g g + gi g )W p = {D-p + j)W p (C.7) 

which means that {D — p + j)~ l ig is the homotopy operator corresponding to g. From 
(C.7) it follows that if Hp satisfies gHp = 0, then (with the exception of the case in which 
j = and p = D) 

Ul = g((D-p + j)-\U p ) (C.8) 
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The operations d and g satisfy the following relations 



d 2 = 0, g 2 = 0, gd + dg = (C.9) 

which means that they define a (commutative) bi-complex. In the following we will use 

dH = ^=> H = dB (CIO) 
g U = ^=> H = gB (C.ll) 

The second relation has been motivated a few lines above. The first is the already mentioned 
local Poincare lemma. Its on-shell validity relies on the construction of a suitable homotopy 
operator which will be discussed below. 

Let us now apply this formalism to prove a generalization of (C.l). Let us start from 
a p-form 

k 
3=0 

and introduce a more compact notation for (C.12) 

k 

J ? = £ W A(J) (°- 13 ) 

j=0 

which enables us to get rid of all tensor indices in the formalism. In formulas such as (C.13) 
it is understood that are p-tensors. We assume now that Jg is closed, as in (C.3). 
This enforces the following descent equations 

dA<°> = 

dA {j) = -gA^-V , j = l,...,k (C.14) 
= gA<® 

The last equality implies that has to be g-exact, i.e., that there exists a p~J-tensor X 
such that 

A (k) = gX (C.15) 

We want to find a (p— l)-form satisfying (C.l). Using linearity and smoothness in 
£, we can conclude that can be written as 

oo 

Q 5 = E W p(i) ( c - 16 ) 

3=0 
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where PW are -tensors that we want to determine by using (C.3). Inserting (C.16) and 
(C.13) in (C.l) we obtain 

dp(°) = A(°) (C.17) 
gP {j ~ 1] + dP {j) = , j = l,...,k-l (C.18) 
ffPM+dP^ = gX (C.19) 
gP ij ~ 1] + dP ij) = , j = k + 1, . . . (C.20) 

In (C.19) we have used (C.15). Observe that from (C.19) it follows that P( fc ) must have 
form 

pW = dZ + gZ' (C.21) 
Now we use the fact that solution to the descent equations is not unique due to the freedom 

Q 5 ~ Q f + dR 5 (C.22) 
where R^ is arbitrary (p — 2)-form. We take 

oo 
3=0 

where are p_ 2 -tensors we can freely choose. Therefore 

oo 

dR ? = ^2(d j (sR " 1 ) + dR (i) ) , (C.24) 

3=0 

which means that we can freely redefine so that 

pQ) „ P (i) + ^rO- 1 ) + dR(j) (C.25) 

It then follows from (C.25), (C.21) and (C.20) that by properly fixing R^', j > 1, we can 
make 

p0') = , j = k,k + l,... (C.26) 

This means that in (C.16) one can terminate the sum with j < (k— 1), as claimed by Wald. 

To avoid unnecessary complications, we will explicitly solve the descent equations 
(C.17)-(C.20) in the case we are primarily interested in, i.e., when p = D — 1 and k = 2 in 
(C.13). Then the equations become 







A (0) 


(C.27) 


<?p (0) 


+ dpV = 


A (i) 


(C.28) 




<?P (1) = 


sx 


(C.29) 


From (C.29) it directly follows that 








P (i) 


= X + g(. 


■•) 


(C.30) 
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But from (C.25) it is easy to see that by appropriately choosing R,(°) we can discard g-exact 
terms in (C.30), so we are left with 

P (1) =X (C.31) 

Now we use (C.31) in (C.28) to obtain 

gP (0) = A (D _ ^ 

which we further operate with i g from the left. Using (C.7) on the LHS we finally obtain 

P(°) = \ig{& {X) -dX) (C.32) 

(C.31)-(C.32) together constitute the searched for solution to (C.l). 

Let us pause here to check consistency of this solution. This means that we must 
obtain (C.27)-(C.29) from (C.31)-(C.32) using only the dJ = descent equations(C14). 
Relation (C.29) is a trivial consequence of (C.31). Then, to obtain (C.28) we operate on 
(C.32) from the left, which gives 

ff p( ) = i^(A«-dX) 

= A (1) - dX - i g g{A- {1) - dX) 
= A (1) - dX - i g {gA {1) + dgX) 
= A« - dX - igigA® + dA^) 

= A« - dpW (C.33) 

and the last line is exactly (C.28). (In the second line we used (C.7), in the fourth line we 
have used the j = 2 relation from (C.14), and in the fifth line (C.31)). Finally, to obtain 
(C.27) we operate with the exterior derivative d on Eq. (C.33) from the left, which gives 

dgPW = dA^ = -gA^ (C.34) 

In the last equality we have used the j = 1 relation from (C.14). Using dg = —gd and 
(C.ll) we obtain 

dp(°) = A^ + g(. . .) = A<°> (C.35) 

which is exactly (C.27). (In last equality we used the fact that (. . .) would have to be 
^_ 2 -tensor, which does not exist (it is by definition identically zero).) This completes the 
consistency check. 

C.2 Contribution to from the pure gravitational C-S Lagrangian term 

We now apply the formal solution of the previous subsection to the contribution of pure 
gravitational CS Lagrangian term to J^. We have 

J ? = 0COV + nc _ i?LpCs _ H? 
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where one uses (1.4), (2.6), (2.14), and (3.10). Obviously Jg has the form (C.12) with k = 2 
and p = D — 1 = 2n — 2, so from the previous subsection we know that we can take to 
have the following form 

Q c = («9£)P (1) + £P (0) = Q$ ) + Qf ) (C.36) 

where P^) are determined by (C.31)-(C.32). Let us calculate them explicitly. 

From 0j? c (i5r) — we can extract X, which when used in (C.31) and (C.36) gives the 
part of proportional to A a b = db^ a 

I' 1 

Qf'=n(n-1) dtP n (A,r,R™- 2 ) (C.37) 
Jo 

Though the calculation of is in principle equally straightforward, it is technically 
much more involved due to the more complicated formulae (C.32) and also because there 
are several terms in (C.36) which contribute to A^ 1 ). The final result is 

(QfV-A^-3 =n(n-l) f dtPnii^TuB^- 2 )^...^ 

Jo 

-n& Jf°V.. WB _ s + ^(2n - iK^t^-a] (C.38) 

C.3 Contribution to of pure gravitational C-S term (alternative derivation) 

Since we know that the term linear in £ (that is Ql ) of is irrelevant for the purpose 
of calculating the entropy, we can follow a simplified path. Starting from 

J 5 = 0- + 0-_^L p cs-H € 

we have already proved in Appendix B, eq.(B.9,B.8), that can be written 

J ? = 0f v _ nP n {i^T + A, R n_1 ) + dqg (C.39) 

where 



q^(r) = q u (r) + q 2 , c (r) 



n(n — 1) 



f 1 dtP n (A,T,-R?- 2 )+ f dtP n ^T,T 
Jo Jo 



(C.40) 



Let us concentrate on 0^ ov — nP n (i^T + A, R n 1 ). We know that this form must be closed 
on shell, and it is linear in £ and dt; only. Therefore we can write it in the form 

Qf v - nP n (i^T + A, R"" 1 ) = £A W + (C.41) 

It follows that 

dA<® « 

gA® + dA« « 

SA« « 



- 26 - 



We expect of course 

ef v - nP n ^T + A, R"" 1 ) « ci(£P (0) ) (C.42) 

because, based on the argument corresponding to eq.(C.26), we can set P( J ') for j > 1. 
Moreover 

P<°) » i^A^ (C.43) 

So (C.39) can be written simply 

J € ps dQ ? (C.44) 

with = + £,i g A^ and given by (C.40). For our later purposes we do not need to 
compute A^ 1 ) explicitly. 

C.4 The variational tri-complex 

In order to justify the use of the d, g complex, introduced in section 2, on shell, i.e. with 
the specification 'modulo the equations of motion', the simplest way is probably to embed 
such complex in the so-called variational complex. This complex is defined by means of 
the operators dn and dy. Consider a space spanned locally by coordinates x , . . . , x n and 
a set of fields and their derivatives 

Aa) = , («) _ aV (Q) 

vi 9W..i fc dx^dx^ ...ax**' 

and define the generic functional derivative 

. . . U l/o' /J 

# = AT ~* = fc , %,/ (C45) 

^(ili2— ifc) 

where /j is the number of occurrences of the integer j amongst ii«2 ■ ■ - ik- Define moreover 
the contact 1-forms 



) - <t>f)dx> (C.46) 



Then for any function or form / made out of a set of fields and their derivatives, we 
have 

dnf = (J^ + Paf)^ + + • • •) dx* (C.47) 

and 

k 

<*V= (C48) 

|7|=0 

du is the form the de Rham exterior operator takes in the context of local filed theory, dy 
is the mathematically consistent way to represent field variations. 
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Then we have 

d% = 0, d\ = 0, d H d v + d v d H = (C.49) 

These two operators define a bi-complex by acting on the space of forms fip of order p in 
space and of order r in the vertical direction (that is, proportional to the exterior product 
of r field variations) . This is known as the variational bi-complex, [29] . 
Now, it is easy to show that 

dug + gdu = 0, d v g + gd v = (C.50) 

thanks to the fact that the operation g corresponds to inserting a constant Kronecker delta. 

We can enlarge the bicomplex to include g and form a commutative variational tri- complex. 

[r j) 

dn,dy and g act on the spaces of forms O v , where j is the same as in section 2. Here 
we reproduce only the relevant (for us) corner of the complex 



n 



(1,2) 




n 



(0,2) 
n-2 



n 



(1,2) 
n-1 



'hi 




n 



(0,2) 
n-1 



" in 




O (0'2) 



■jr(l,2) 




/ 9 / / 9 




The map J is the projection that maps out total differentials. The above diagram is 
commutative and the complex is (locally) exact. The exactness for dn and dy for the bi- 
complex is proved by constructing the corresponding homotopy operators. The homotopy 
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operator for g has been constructed in section 2, (C.7). Therefore the tri-complex is exact. 
The spaces T > 3 are spaces of equations of motion. In fact the space of equations of 
motion splits into subspaces labeled by j. For instance, 6^g = £,g^ + d^g^. Therefore 
5^g E = £E( ) + <9£E ( - 1 ). But E^ ) and E*- 1 ) are proportional, so they vanish simultaneously 
on shell. This entitles us to split the tri-complex according to j and treat the corresponding 
pieces separately, as we have done in subsection C.l. 



D Geometry of codimension-2 surfaces 

We review here some basics of the geometric formalism of codimension-2 surfaces ((D — 2)- 
dimensional surfaces embedded in D-dimensional space) which we apply to a special case 
of bifurcation surface in Section 4.2. We follow the approach reviewed in [30, 31], and 
especially Section 2 of [32]. 

On a codimension-p surface S embedded in a D-dimensional Lorentzian spacetime M 
one can separate the metric g a b of Ai in normal and tangent parts 

9ab = Kb + q a b (D.l) 

where q a b is known as the first fundamental tensor or induced metric on S. In our ap- 
plications S is spacelike, so the tangential part q a b is Riemannian, and the normal part 
Kb is Lorentzian. A decomposition (D.l) is obtained by demanding for q a b and h a b to be 
projection operators satisfying 

q a b t b = t a , q a b s b = 0, h a b t b = 0, h a b s b = s a (D.2) 

for an arbitrary vector t a tangential to 5, and an arbitrary vector s a normal to S. It follows 
that q a b and h a b satisfy the standard projector relations 

q a b q b c = q a c h\h a c = h a c , q a b h b c = (D.3) 

If £ a and n a is an arbitrary pair of two (future directed) null vector fields satisfying t a Ti a — 
— 1, then we can write the normal part of the metric as 

Kb = -£ a n b - n a £b (D.4) 

One can express h a b in terms of the binormal e a b of S. As binormal can be written as 

e a b = t a n>b ~ nJb (D-5) 
From (D.4) and definition of n a and £ a it follows that 

Kb = ta^cb (D.6) 

In our application in which S = B we can take n a = v a and £ a = u a , where u a and v a are 
vector fields obtained from Kruskal-type coordinates and defined in (F.7) in Appendix F. 
Note that generally l a and n a are not uniquely defined, and the freedom is in the rescaling 

T -»• ft , n a -»• n a /f (D.7) 
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where / = f{x) is an arbitrary smooth positive function. 

Another important object is the second fundamental tensor K ab c which can be ex- 
pressed as 

K ab c = q a d q b V d q e c = -q a d q b e V d h e c (D.8) 

where the second equality follows from (D.l) and V a 9bc = 0. 

Of our main interest is the normal curvature tensor Q ab which can be obtained from 

n ab = ~\e cd (q e a qlh9 c h h d R efgh + K aed K b e c - K bed K a e ^ 

= ~ \<f a <lU cd Refcd ~ \e cd (K aed K b e c - K bed K a e c ) (D.9) 

There is a connection u) a , the connection of the 50(1,1) normal bundle, associated with 
the normal curvature tensor defined by 

u a = -q e a n d V e £ d (D.10) 

It is important to emphasize that uj, being a connection, is pseudo 1-form, which means 
that it depends on the choice of the "frame" vectors l a and n a . It is easy to check that 
under the rescaling (D.7) (which is related to change of frame) it transforms as 

co a ^LO a + q b a V b lnf (D.ll) 

which can be viewed as a gauge transformation of an £0(1, 1) connection. Note that h ab , 
q ab , K ab c and VL ab are instead normal (not pseudo) tensors, i.e., they are invariant on the 
gauge transformation (D.ll) (because they all can be defined without using frame vectors 
l a and n a ). 

It can be shown that for codimension-2 surfaces one gets 

^ab = q a qii^cUd ~ V d w c ) = q c a qt{duj) cd (D.12) 

In the case of our interest, when S is bifurcation surface B, from (D.l) and (D.6), due 

B 

to the special property V a e a {, = 0, it follows that 

V a h bc = , V a q bc = (D.13) 
Using this in (D.8) implies that the second fundamental tensor vanishes on B 

K ab c = (D.14) 

Using this in (D.9) we finally obtain 

Slab = -\q e a qU cd Refcd = q e a q f b (RN)ef (D.15) 
where the 2-form R^r was defined in (4.18). 
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E Entropy formula and spin connection 

The purpose of this section is to outline a derivation of a covariant formula for the entropy 
in a coordinate independent way. We defer a proof, which needs the usage of the vielbein 
formalism, to a future paper (for a Lorentz invariant formulation in 3D, see ref.[6]). In order 
to appreciate the following reduction procedure it is useful to review the geometrical setting 
imposed upon us by the problem we are studying, see [33], vol.11. The geometry is that of an 
asymptotically Minkowski space time manifold X with a codimension 2 submanifold B. We 
have O(X), the bundle of orthonormal frames on X with structure group SO(D — 1,1) and 
0(B) the bundle of orthonormal frames on B with structure group SO(D — 2). We consider 
also the bundle of adapted frames. An adapted frame is a complete set of orthonormal 
vectors which are either tangent or orthogonal to B. They form a principal bundle 0(X, B) 
with structure group SO(l, 1) x SO(D — 2). To complete the description we have the bundle 
of normal frames ON(B) with structure group 50(1, 1) and the embedding i: 0(X, B) 
0(X). For convenience, let us denote by f) and 5 the Lie algebras of SO(D — 2) and 
50(1,1), respectively. 

Let us come now to formula (4.11). The connection T in it is a connection of the linear 
frame bundle LX. Every metric connection in LX is in one-to-one correspondence with a 
connection in 0(X) (see [33], vol.1, ch. 4, § 2). So we assume we can replace (4.11) with 

5 CS = 2?m(n - 1) C dt [ P n (e, a, R™~ 2 ) (E.l) 
io Jb 

where a is the reduction of T to the structure group 50(1, D — 1) and R its curvature. 
Moreover e = EeE -1 , where E = {E 1 ^} are the vielbein in X (a is a flat index, see below). 
P n is now a symmetric polynomial in the Lie algebra of SO(D — 1,1). 

In (E.l) it is understood that the forms in the integrand are pulled back from X to B. 
Now by pulling back a generic connection a of 0(X) through i, we do not get a connection, 
unless we restrict to the components in f)+ t. If so, the connection splits into at + a m that 
is a connection at in 0(B) with values in f) and a connection a n in ON(B) with values in 
t (see [33], vol.11, ch. VII). As we have clarified in section 4, the geometry of the problem 
is defined by the presence of the surface B with its tangent and normal directions, thus the 
just considered reduction of a connection pulled back from X, is natural in this scheme. 
But once we replace in (E.l) the connection at + a n , with values in the direct sum f)+ fi, 
the presence of the binormal e maps out the F) components and only the components along 
t (the Lie algebra of the normal frame bundle structure group 50(1, 1)) survive. At this 
point we have to do with an Abelian connection and we can easily integrate over t as in 
section 4.1. 

To view the situation in more detail let us introduce the following conventions. Let us 
denote by fi, v, .. = 0, . . . , D — 1 generic curved indices and a,b, . . . = 0, . . . , D — 1 generic 
flat indices. Then, following [30], we will denote by A, B, .. = 2, . . . ,D — 1 flat tangent 
indices in B and by X, Y, . . . = 0, 1 normal flat indices, and introduce adapted vielbein ia^ 
and Ax M (they are particular cases of E%) , such that 

<? v = ia»i A v , W v = \x fl X X u (E.2) 
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One can show in particular that, since Ax^Ay^ = r\xy (tj denotes the flat Minkowski 
metric), one can make the following identifications 

V-^f. V-^ (E.3) 

with reference to the null vectors introduced in the previous Appendix. Then, it is easy to 
show that 

ep,E£E% = rj la rio b - r] lb r] 0a (E.4) 

where r/ is the flat Minkowski metric. Thus, for instance, 

tr(ea n ) = e^E p E^(a n ) ab = 2af (E.5) 

and likewise for the curvature. Therefore in this approach we obtain formulas like (4.18) 
and (4.17) with T and R replaced by a n and its curvature R n . It is understood that all 
the forms are pulled back to B, which can be achieved on components by contracting the 
form index with the q projector: for instance the intrinsic component of the pulled back 
a n is 

It is now convenient to compare the normal bundle connection with the one introduced 
in [30], 

mjp = K v \ X p V li \x (T , where V fl = q%V u (E.6) 
Using VEa = —a a b E^ we can rewrite 

«Vp = g/e/C&X (E.7) 

Saturating with e u p we obtain precisely the RHS of (E.5). 
On the other hand, inserting (E.3) into (E.6) one finds 

ro/p = -e> r V M r (E.8) 

Saturating with e p u and dividing by 2, we get precisely the definition (D.10). 

Finally a comment about gauge transformations in the normal frame bundle. They 
are valued in £0(1, 1) and act on Iq, l\ as follows 



Iq \ I cosh t sinh t\ ( Iq 
lij \ sinh t cosh tj \l\ 



(E.9) 



where t is a local function. Using again (E.3), it is easy to see that they act on n,l as a 
rescaling 

n -> e*n, I -> e~H (E.10) 
Identifying 1/ f = e we get precisely the rescalings considered in (D.ll). 
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F Kruskal-type coordinates 



Here we prove some of the relations and properties used in Section 4. The strategy we 
use is to first make calculations in special "Kruskal-type" coordinates, and then eventually 
generalizing to other coordinate systems typically used in black hole calculations. 

In [26] it was shown that, in a spacetime with Killing horizon on which the surface 
gravity is constant, one can construct Kruskal-type coordinates (U, V, {x 1 }), i = 1, . . . , D — 
2, in which metric has the following form 

ds 2 = G dUdV + VHi dx l dU + g {j dx i dx j (F.l) 

where G, Hi and gij are generally smooth functions of D — 1 variables U, V and {x 1 }. The 
physical horizon is at U = 0, while U = V = defines the bifurcation surface B. On the 
bifurcation surface B we have G\b = —2/k where k is the surface gravity and is constant 
throughout B. We see that {x 1 } are tangential and U,V are normal on B. The horizon 
generating Killing vector field £ is given by 

£ = K (uJL-V-^\ (F.2) 

where the constant K is surface gravity. In Kruskal coordinates, the components of the 
metric (F.l) are regular and well-defined on B, and the components of £ obviously satisfy 

r| B = o , v b e\ B = d b e = A-\ (f.3) 

and the nonvanishing components of A a b are 

= -A v v = k (FA) 

From this it follows that 

dA a b = d M AV^ = (F.5) 

Using Kruskal-type coordinates one can introduce a pair of null vector fields and 
normal to B defined by 

"-Or)'- 



If we make a rescaling 



(F.7) 



then and are a pair of future-oriented null vector fields normal to B which is 
normalized such that 

Ufl v» = -1 (F.8) 



- 33 - 



Using Kruskal coordinates we can obtain important information on the components of 
connection and Riemann tensor on the bifurcation surface B. We start with the Killing 
lemma 

V a V 6 f = R c ba dZ d (F-9) 
satisfied by any Killing vector field. On B the left hand side of (F.9) can be written as 

V a V b e 1 T c hp A? a - r ab A\ + T c as A\ (F.10) 

so that, 

n p A? a -T r ab A c r + T c as A s b 1 (F.ll) 
By making particular choices for free indices in (F.ll) we get 

F^ = Tfj = Ty i = T^j = Vju = T l j V = T\ju = Ty V = (F.12) 

and 

Tuu = ^vv = Tyy = ^uu = Tuv = ^uv = (F.13) 

B 

All this equalities should be understood as =, which, for brevity, we avoid to write here. 
It follows that the only nonvanishing components of the pull-back of the tensor valued 

Uii 1 Vi ana 1 ji- 



1-form connection T on the bifurcation surface B are T^-, Ty- and T^-. One important 



consequence is 

T a ci T c bj dx i A dx j = , when a and/or b is U or V . (F.14) 
By taking the covariant derivative of (F.9) we get 

V p V a V b r 1 R c bad A d p = (r c bdta - Ti a4 + T ae Tt d - T df TQ A d p (F.15) 
while from the definition of covariant derivative we have 

%v a v fe e c = rg dja A d p + r a£jP a\ + rg /iP Af a - r^ p A c g + 

+ (T ar Tl h -r s ab r c sh )A h p (F.16) 
Equating right hand sides gives us 

n a ,d K + Ke, P A e b + n f , P A J a - Kb, P A C 9 = (F.17) 
Furthermore, from (F.17) we have 

^Yj,k = ^Yj,k = ^Vi,j = ^Ui,j = ^Uj,k = ^Vj,k = (F.18) 

B 

Again, all the above equations should be understood as = . 
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Now we focus on the Riemann tensor R a b ^ lll - As we do the integration over the bi- 
furcation surface, the only relevant components are those for which the "form indices" ^ 
and v relate to coordinates from the subset {x 1 }. Using results from above inside the the 
definition of Riemann tensor, 

rya -pa _ -pa , -pa ptr -pa -pp /-p i q\ 

-n- bpv ~ 1 bu,p 1 bp,u ' L pa 1 bv 1 vp L bp \ s > 

we obtain 

R U kij = , R k uij = > R V kij = i R k Vij = (F.20) 

In other words, the 2-form R a bi j is block-diagonal in the tensor indices when pulled-back 
on the bifurcation surface B. 

We can write the obtained properties of V&£ a and Riemann tensor in the following 
covariant way 

= q c b V c C 1 q a c V b e 1 q c b e\ 1 q\e c b (F.21) 

1 h\ qh q\ q r d R e fgr 1 q\ h? b q 9 c q r d R e fgr (F.22) 

where q a b and h a b are tensors obtained by separation of metric tensor g a b into tangent and 
normal part to B, defined in Appendix D. Being written in a covariant way, relations in 
(F.21) and (F.22) are valid in all coordinate systems. 



G Schwarzschild-type coordinates for stationary black holes 

Here we give the proof of (4.30) in the case of stationary geometry. The relation between 
Kruskal-type coordinates (U, V, {x 1 }) (described in Appendix F) and Schwarzschild-type 
coordinates (i,r, {x 1 '}) was constructed in [26], and is given by 



t = — In 

2k 



-^exp {2KH{UV,{x i })) 



V 

r = UVf(UV,{x 1 }) (G.l) 
x l = x l , i = 2, . . . , D — 2 (coordinates tangent to B) 

where / and H are some smooth functions around B (defined by U = V = 0), which are 
non- vanishing on B. 

The difference between Tn evaluated in two (arbitrary) coordinate systems {x a } and 
{x a '} is 

Ar> = ^triedAA" 1 ) = -itr(eAcZA^) (G.2) 

where (A) a b = d v x a and (A" 1 ^ = d b x a ' . 

When we specialize to {x a } = (U, V, {x 1 }) (identified below by index K, meaning 
Kruskal-type) we obtain 

2k K T N = (A K dA K ) v v - (A K dA K ) U u (G.3) 
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There is no reason for the expression in (G.3) to be vanishing in general. However, we shall 
now show that this happens when we take {x a '} = (t, r, {x*}), i.e., Schwarzschild-type. First 
of all, note that integration over B forces form-indices to take values i, i = 2, . . . , D — 1, so 
we need to calculate 

dV d 2 x a ' dU d 2 x a ' 
2A K (T N ) l = dxa , dvQx . - Qxa , duQx . (G.4) 

From (G.l) it follows that the nonvanishing second derivatives present in (G.4) are^ 
d 2 t ■ d 2 t 



dU dx l dVdx 1 

wkf - v ( Btf + UVBJ ) ■ wk - u + UVBtf 

where a dot on a function means g = duv9- We also need 

dU dU 1 dV dV 1 

where /(0) = f(0,{x t }). By using all this in (G.4) we obtain that the lowest order terms 
in (G.4) cancel and 

A K (T N )i = 0(UV) A (G.6) 

which implies (4.30). 

Note that in coordinates (t,r,{x 1 '}), the horizon generating Killing vector £ is equal 
to d/dt. This means that (G.l) does not include coordinates of the Boyer-Lindquist type 
where £ is of the form £ = d/dt + Qd/d(j). Fortunately, it is straightforward to extend 
the result A^(rV)i = to the following more general transformations which include 
coordinates of the Boyer-Lindquist type: 



t = — In 
2k 



-^-exp(2 K H(UV,{x i })) 



V 

r = UVfXUV^x 1 }) 
1 = x l (U, V, {x- 7 }) , i,j = 2,...,D — l (x J are coordinates tangent to B) 



where the Jacobian 



References 



dx a 



dx b 



7^ is well defined on B (a, b = 0, . . . , D — 1). 



[1] S. Deser, R. Jackiw and S. Tcmplcton, Three-dimensional massive gauge theories, 
Phys. Rev. Lett. 48, 975 (1982) 

[2] S. Dcscr, R. Jackiw and S. Templeton, Topologically massive gauge theories, Ann. Phys., NY 
140, 372 (1982) 



**As already commented before, when calculating in singular coordinate systems we have to avoid using 
U — V = in the middle of calculations. This is why we are calculating for finite U, V, and only at the end 
of the calculation we take U, V — > 0. 



- 36 - 



[3] E. Witten, (2+1) Dimensional Gravity as an Exactly Soluble System Nucl.Phys. B311 
(1988) 46. 

[4] S. Deser and B. Tekin, Energy in topologically massive gravity, Class. Quantum Grav. 20, 
L259 (2003) [arXiv:gr-qc/0307073vl] 

[5] A. Maloney and E. Witten, Quantum Gravity Partition Functions in Three Dimensions, 
JHEP 1002, 029 (2010) [arXiv:0712.0155 [hcp-th]]. 

[6] S. N. Solodukhin, Holography with gravitational Chern-Simons, Phys. Rev. D74 (2006) 
024015. [hep-th/0509148]. 

[7] R. F. Perez, Conserved current for the Cotton tensor, black hole entropy and equivariant 
Pontry agin forms, Class. Quant. Grav. 27 (2010) 135015. [arXiv:1004.3161 [gr-qc]]. 

[8] P. Kraus, F. Larsen, Holographic gravitational anomalies, JHEP 0601 (2006) 022. 
[hep-th/0508218]. 

[9] M. -I. Park, BTZ black hole with gravitational Chern-Simons: Thermodynamics and 
statistical entropy, Phys. Rev. D77 (2008) 026011. [hep-tli/0608165]. 

[10] O. Miskovic, R. Olea, Background- independent charges in Topologically Massive Gravity, 
JHEP 0912 (2009) 046. [arXiv:0909.2275 [hep-th]]. 

[11] S.S. Chern, Complex manifolds without potential theory, Berlin, Heidelberg, New York: 
Springer 2nd cd. 1979. 

[12] A. Sen, Black Hole Entropy Function, Attractors and Precision Counting of Microstates, 
Gen. Rcl. Grav. 40 (2008) 2249 [arXiv:0708.1270 [hcp-th]]. 

[13] P. Dominis Prcster and T. Terzic, a' -exact entropies for BPS and non-BPS extremal dyonic 
black holes in heterotic string theory from ten- dimensional super symmetry, JHEP 0812 
(2008) 088 [arXiv:0809.4954 [hep-th]]. 

[14] P. Dominis Prester, AdS^ backgrounds from 10D effective action of heterotic string theory, 
Phys. Rev. D 81 (2010) 044014 [arXiv:0912.0030 [hcp-th]]. 

[15] S.S. Chern and J. Simons, Characteristic forms and geometric invariants, Ann. Math. 99 (1) 
(1974), 48. 

[16] L. Bonora, P. Cotta-Ramusino, M. Rinaldi and J. Stashcff, The evaluation map in field 
theory, sigma-models and strings. I, Comm. Math. Phys. 112 (1987) 237. 

[17] V. Iyer and R.M. Wald, Some properties of the Noether charge and a proposal for dynamical 
black hole entropy, Phys. Rev. D 50, 846 (1994) [arXiv : gr-qc/9403028vl] 

[18] Y. Tachikawa, Black hole entropy in the presence of Chern-Simons terms, Class. Quantum 
Grav. 24, 737 (2007) [arXiv:hep-th/0611141v2] 

[19] S. N. Solodukhin, Holographic description of gravitational anomalies, JHEP 0607 (2006) 003 
[arXiv:hep-th/0512216]. 

[20] A. A. Garcia, F.W. Hehl, C. Heinicke and A. Macias, The Cotton tensor in Riemannian 
spacetimes, Class. Quantum Grav. 21, 1099 (2004) [arXiv : gr-qc/0309008v2] 

[21] R.M. Wald, Black hole entropy is the Noether charge, Phys. Rev. D 48, R3427 (1993) 
[arXiv : gr-qc/9307038vl] 

[22] T. Jacobson, G. Kang and R.C. Myers, On black hole entropy, Phys. Rev. D 49, 6587 (1994) 
[arXiv : gr-qc/9312023v2] 



-37- 



[23] R.M. Wald, On identically closed forms locally constructed from a field, Jour. Math. Phys. 31 
(1990) 2378. 

[24] R. M. Wald, A. Zoupas, A General definition of 'conserved quantities ' in general relativity 
and other theories of gravity, Phys. Rev. D61 (2000) 084027. [gr-qc/9911095]. 

[25] I. Racz and R. M. Wald, Global extensions of space-times describing asymptotic final states 
of black holes, Class. Quant. Grav. 13 (1996) 539 [arXiv:gr-qc/9507055]. 

[26] I. Racz and R. M. Wald, Extension of space-times with Killing horizon, Class. Quant. Grav. 
9 (1992) 2643. 

[27] L.Bonora and M.Rinaldi Normal bundles, pfaffians and anomalies, Nucl Phys. B 578 (2000) 
497 

[28] M. Nakahara, Geometry, Topology and Physics, Taylor & Francis (2003) 
[29] I.M.Anderson, The variational bicomplex and references therein. 

[30] B. Carter, Essentials of classical brane dynamics, Int. J. Theor. Phys. 40 (2001) 2099 
[arXiv:gr-qc/0012036]. 

[31] I. Booth and S. Fairhurst, Isolated, slowly evolving, and dynamical trapping horizons: 
geometry and mechanics from surface deformations, Phys. Rev. D 75 (2007) 084019 
[arXiv:gr-qc/0610032]. 

[32] L. M. Cao, Deformation of Codimension-2 Surface and Horizon Thermodynamics, 
arXiv:1009.4540 [gr-qc]. 

[33] S.Kobayashi and K.Nomizu, Foundation of differential geometry , vol.1, II, Johm Wiley & 
Sons, New York 1963. 

[34] D.D.K. Chow, C.N. Pope and E. Sezgin, Classification of solutions in topologically massive 
gravity, Class. Quantum Grav. 27, 105001 (2010) [arXiv : 0906 . 3559v3] 



- 38 - 



